On the relative positions of the 2A peaks in Raman and tunneUng spectra of d— wave 

superconductors. 
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We study Big Raman intensity i?(Sl) and tiie density of states N{uj) in isotropic 2D d— wave 
superconductors. For an ideal gas, R{Q) and N{uj) have sharp peaks at f2 = 2A and lo = A, 
respectively, where A is the maximum value of the gap. We study how the peak positions are 
affected by the fermionic damping due to impurity scattering. We show that while the damping 
generally shifts the peak positions to larger frequencies, the peak in R{Q,) still occurs at almost twice 
the peak position in Af(w) and therefore cannot account for the experimentally observed downturn 
shift of the peak frequency in R{Q,) in underdoped cuprates compared to twice that in N{u)). We 
also discuss how the fermionic damping affects the dynamical spin susceptibility. 



PACS numbers;71.10.Ca,74.20.Fg,74.25.-q 

The unusual physical properties of cuprate supercon- 
ductors have continued to be of high interest to con- 
densed matter physicists for more than a decade. In re- 
cent years, a lot of attention was devoted to the study of 
collective bosonic excitations in the superconducting and 
pseudogap phases ||l|-^ . The most notable experimental 
observation is the discovery of the resonance mode in the 
dynamical spin susceptibility QJ^. This mode is cen- 
tered at the antiferromagnetic momentum Q — (tTjTt), 
and physically reflects the fact that near antiferromag- 
netic instability, collective spin excitations in a d-wave 
superconductor are undamped, propagating spin waves 
at energies smaller than 2 A a-M- Less attention is de- 
voted to the study of possible resonance bosonic excita- 
tions at zero momentum transfer. These excitations are 
probed by Raman scattering which generally measures 
the imaginary part of the fully renormalized particle-hole 
susceptibility at vanishingly small incoming momentum, 
weighted with Raman form factors which depend on the 
scattering geometry pO| , pd| . The experiments relevant 
to our discussion were performed in Big geometry where 
the Raman form factors are the largest for fermionic mo- 
menta near (0, tt) and symmetry related points where the 



d^2_y2 gap A(A;) is near its maximum A |12|Jl 

In a BCS theory for a d— wave superconductor Big Ra- 
man intensity i?(J7) logarithmically diverges at 2A and 
rapidly, as to^, decreases at smaller frequencies [|3|. Ex- 
perimental data for overdoped Bi22l2 are qualitatively 
consistent with this behavior 11,0. Furthermore, the 2 A 
extracted from _R(57) is almost exactly twice the gap ex- 
tracted from SIN tunneling data, which measure a single 
particle density of states (DOS) N{lj) 0. With under- 
doping, however, the peak frequency in R{fl) progres- 
sively deviates down from the 2A extracted from the 
tunneling experiments 1^. 

Blumberg, Morr and one of us (CBM) |15| attributed 
this deviation to a final state interaction between scat- 
tered quasiparticles. They argued that the magnetically 
mediated final state interaction in Big geometry is at- 



tractive and gives rise to a pseudo-resonance in R{n) at 
a frequency ^res, which with underdoping progressively 
deviates down from 2A. 

An alternative to the resonance mode scenario is one 
in which the final state interaction is irrelevant, and the 
discrepancy between Raman and tunneling data is due 
to fermionic incoherence, which generally shifts the posi- 
tions of both the Raman peak and the peak in the DOS. 
In this paper, we show that the shifts in the peak posi- 
tions of i?(r2) and N{u) due to fermionic damping are 
correlated such that without final state interaction, the 
peak in R{^) is still located at almost exactly twice the 
peak frequency in N{iu). This result implies that the 
experimentally observed relative downturn deviation of 
the peak in the Raman intensity cannot be explained by 
purely fermionic self-energy effects, and leaves the reso- 
nance mode scenario as the most probable one. 

We begin with the general expressions for R{fl) and 
N{uj) in a superconductor. The DOS is the imaginary 
part of the local normal quasiparticle Green's function, 
and the Raman intensity without final state interaction 
is the imaginary part of the fermionic polarization bubble 
with a zero momentum transfer, weighted with Raman 
vertices pO,Gl[. We have, up to an overall factor 



R{n) = Im / dk duj 



Vb 



(k) {Gscik,uj+) Gsc{k,uj-) + 



F{k,uj+) F(fc,tj_)); N{uj) ^Im dk Gsc{k,oj) (!) 

Here, Vbi (k) ex cos k^ — cos ky is the Big Raman vertex, 
u;± = u;±n/2, and Gsc{k, lj) and F(k, lu) are normal and 
anomalous quasiparticle Green's functions given by 

G^k, oj) - G-\-k, -Lo)/{G-\k, u) G~\-k, -L,) + Al) 
Fik,u;)^iAk/iG-\k,u;)G-\-k,^cj) + Al) . (2) 

Here, G~^{k,uj) — I](k, tj) — e^, where I](k, w) is the 
fermionic self-energy (which also absorbs a bare lo term) , 
and Afc is the superconducting gap, which for a dx2_y2 
superconductor behaves as A^ oc coskx — cos ky. 



In an ideal gas, the fermionic self-energy is absent (i.e., 
S(k, w) = uj). To simplify the discussion, we assume 
that the Fermi surface is circular, for which both R{U) 
and N(lli) can be evaluated exactly |1^. Substituting 
the momentum integration by integration over dek, and 
approximating the fc— dependences of Vbj^ (k) and A^ by 
cos 8, we obtain for fi, a; < 1 



R{n) 



^(73^1,^,3,(^2), iV(c^) 



u,K{u^), (3) 



and for Q, u > I 



R{n)^ 



37r ,15 1 , , , 2 /I 



(4) 



where w — uj/A, ft — il/{2A), and F{a, b, c, x) and K{x) 
are hypergeometric and elliptical functions, respectively. 





n/2A,(0/A 
FIG. 1. The behavior of the Raman intensity R{i}) and 

the DOS N{Li}) in d— wave superconductors, a), b) Fermi-gas 

results (7 — 0) SoUd lines -d— wave results, dashed lines 

-s— wave results shown for comparison. The insert shows the 

behavior of R{Q,) close to the threshold frequency 2A. c) The 

results for 7 — 0.2A (Eq. (M)) The normalized frequencies 

are ui = uj/A and Q = f2/(2A). Observe that for d— wave 

superconductors, the peak in the Raman intensity is located 

at larger normalized frequency than the peak in the density 

of states. 

The results for R{il) and N{uj) are plotted in Fig. 
nk,b. At the smallest frequencies, R{fl) ex fl^ and 
N{uj) ex cj. At larger frequencies, the Raman inten- 
sity and the DOS diverge logarithmically a,t u; — A 
and n — 2 A, respectively. At even larger frequen- 
cies, the momentum dependence of the gap progres- 
sively becomes less relevant, and both i?(f7) and N{llj) 
acquire the same forms as in s— wave superconductors: 
R(n) ex l/{nVn^ - l), N{u) = w/Vtl-^ _ i. observe 
that R{Q) crosses over to the s— wave behavior immedi- 
ately above 2 A. (see inset in Fig. |k). 



Our goal is to study how the peak positions and, more 
generally, the functional forms of R{fi) and N{uj) are 
affected by the fermionic self-energy. In general, the 
fermionic self-energy comes from various sources, and 
at least part of it, associated with the scattering by 
the same bosonic excitations which give rise to super- 
conductivity, has to be determined fully self-consistently 
from the Ehashberg-type equations [|j. In this paper, 
we assume for simplicity that the primary source for the 
fermionic damping is impurity scattering. We consider 
the fermionic self-energy in the self-consistent T— matrix 
formalism ||l^,|l^ and neglect subtle 2D effects beyond 
T— matrix approximation ||l^. Hirshfeld, Wolfle and 
Einzel |17| demonstrated that for a d— wave superconduc- 
tor with a particle-hole symmetry and fc— independent 
scattering potential, S(fc,a;) = uj + i j^ sign^Lu), where 
7|jj| is a solution of the self-consistent equation |17| , [16[ 



7H = 7 



go(l^l) 

c'+.9o(l^l 



(5) 



Here 7 is proportional to the impurity concentration, c is 
the cotangent of the scattering phase shift, and 50 (1^1) = 
(j/7riV^) ^j, Gsc(fc, (j-") sign(u), where Np is the normal 
state DOS at the Fermi surface. The self-consistency of 
(B) is in the fact that Gsc{k,uj) given by (0) by itself 
depends on 7(0;) through S(w). Physically this means 
that the fermionic self-energy due to impurity scattering 
is by itself affected by a superconductivity. 

For a circular Fermi surface, the momentum inte- 
gral over G{k,Ljj) can be performed exactly and yields 
go{\uj\) = {2/tt)K{A^/Y?{uj)). In the normal state, 
goiu)) = 1, and ^\^\ reduces to a constant ^\^\ = 7/(1 -I- 
c^). In a superconducting state, however, ^\i^\ is complex 
and frequency dependent. Still, one can easily demon- 
strate that 7|^| remains finite for all frequencies. Below 
we will need 70 and 7a- For unitary scattering (c=0), we 
solved (^) at these frequencies and for 7 <C A obtained 
with logarithmical accuracy 7a = (2/7r)7 log(A/7) and 
7o = (7r7A/log(A/7))i/2. 

We now proceed with the calculations of the DOS and 
Raman intensity. As the self-energy is fc-indcpendcnt, 
we can use the same trick as in earlier studies p, p^ , p2[ , 
and first integrate over momenta in (|l|) . Substituting the 
momentum integration by the integration over efc and 
evaluating the integrals, we obtain |l3] 



N{u:) = Im r de ^^'^^ 



^7r/2 



X / du) 
'—00 



Dicj) 

(E+-E_)2 + (Z?+-Z?_)2 



R{n) = -Re / dOcos^e 
Jo 

(6) 



4D+D_(D+ + D_) 

Here uj± = uj±n/2, i:± :^ S(c^±), D± = D{lo±), and 
Diu) = ^A2cos2e-E2(a>). 

As a warm up, consider the limit of small frequencies. 
Substituting S(u;) = uj + ijosign{uj) into (o) and expand- 
ing in frequency, we obtained 



R{n) 

N{u}) 



n 7^ lo, 

N{0)+Q (1 



1/70 + 0(^3) 
2 



-tan 



70^ 

Co' 



270 



logJl + (^)^ 




b) 




^ 



FIG. 2. The behavior of the spin polarization operator in 
the Fermi gas (a) and at a finite 7 (b). Solid line - RenQ{Q,), 
dashed line - ImIlQ{Q,) 

where iV(0) = (2/7r)7o log 1/%, and % = 70/A. At 
(!> < 1, N{u!) = N{0) + <I'V(7r7o). 

We see that fermionic damping (i) yields a linear fre- 
quency dependence of R{fl), and (ii) yields a finite DOS 
at zero frequency, and the quadratic frequency depen- 
dence of N{uj) above N{0). Both of these results agree 
with earlier studies ||l7| , |l6| , |l9|l . 

We now consider R{il) and at fJ « 2A and oj ss A 
where the Raman intensity and the DOS diverge in a 
Fermi gas. Simple estimates show that near 2A, the in- 
tegral in R{i^) is dominated by D± <C 1 i.e. u!± w ±A. 
Substituting E(a; w A) = a; + 17^ sign{Ld) into dq) and 
expanding the integrands to first order in fi — 2A and 
LO — Is. respectively, we obtain after lengthy but straight- 
forward calculations 



R[n) 

N{u)) 



4V2 
V2 



Re 



2' 2 



2, 



r2— l+i7A 



v^-1 



«7A 



Re 



K 



-2 



ti; — l+i7A 



^uj -1 + ijA 



(8) 



Here 7a = 7a/A. 

We now analyse these results. Near ft = Oo = 1, we 
find from (|) R{n) « (-1/8) log[(f^ - 1)^ + 7^^], N{Ld) « 
(—l/47r)log[((D — 1)^-1- 7a]. We see that the logarithmical 
divergencies are cut by the fermionic damping, but, to 
a logarithmical accuracy, the peaks remain at the same 



positions as in a Fermi gas. In other words, fermionic 
damping gives rise to a broadening of the peaks in the 
(' Raman intensity and the DOS, but still, the peak in R{^) 
is located at twice the peak frequency of the DOS. 

Calculations beyond the logarithmical accuracy show 
that the peak positions do shift to higher frequencies but 
the relative shift is opposite to the one detected in the ex- 
periments: the peak in R{Cl) shifts to higher frequencies 
for arbitrary 7a {^peak — 1 oc 7a ^'^^ 7A <C 1), while the 
peak in N^uj) shifts to high frequencies only if the damp- 
ing exceeds a threshold value of 7a « 0.77. Obviously, 
the magnitude of the shift in N{u!) is smaller than that 
in _R(51). The behavior of N{lo) and R{0.) is illustrated 
in Fig. |l|c. 

Note in passing that for s— wave superconductors, the 
the same calculations which lead to (g) yield 



R{x),N{x) 



V^ 



(1 + a2)7A 



(9) 



where ax = {x — 1)/(7a) and x — Cl for R(x) and x — ui 
for N{x). Again, the divergencies at 0. — 00 — I are gone 
and substituted by the peaks at higher frequencies for 
which Ut = —^ ■ At the same time, the functional forms 

of i?(51) and N{Ld) are identical. This implies that in a 
dirty s— wave superconductor, the peak in the Raman in- 
tensity is also located exactly at twice the peak frequency 
in N{uj), though both peak positions shift from the Fermi 
gas values. 

For completeness, we also discuss how fermionic damp- 
ing affects the spin polarization operator at the antifer- 
romagnetic momentum Q. The form of this polarization 
operator is relevant for the interpretation of neutron scat- 
tering and ARPES data |@,i|2|. 

The spin polarization operator is related to the dy- 
namical structure factor as Sq{CI) cx Im{£^~^ ~IlQ{i).))~^ 
where ^ is the magnetic correlation length |g]. It is for- 
mally given by the same set of particle-hole bubbles made 
of normal and anomalous Green's functions as the Raman 
intensity, but differs from R{n) in two aspects. First, 
the antiferromagnetic spin polarization is a finite mo- 
mentum probe, and the contribution to low frequency 
ImIlQ{fl) only comes from the momentum range in the 
Brillouin zone where both particles in the bubble are near 
the Fermi surface hot spots. Near hot spots, the super- 
conducting gap is finite and close to A. In other words, 
the regions near the nodes of the dx2-y2 gap do not con- 
tribute to the dynamical spin susceptibility near Q. Sec- 
ond, the vertices for IIq{U) contain Pauli matrices. For 
the anomalous FF term, the summation over spin pro- 
jections yields an extra factor —1 compared to the Ra- 
man bubble. Performing the momentum integration in 
the GG and FF bubbles in the same way as before, and 
using the fact that AfcAfc+g = — A^, we obtain [g| 



UQin) 



duj- 



E+S- + D+D^ 



2D+D^ 



(10) 



The overall factor is chosen such that in the normal state, 

nQ(0) = i|c^|. 

For an ideal gas, the frequency integration in (|lO| ) 
yields Im UQ{fl) = and Re Uqifl) (x fl^ at fl ^ 
51/ (2 A) < 1. For large enough f , this behavior of IlQ{il) 
gives rise to a resonant peak in Sq{Q) at a frequency 
where ReUgin) = ^^ ^^. At fi = 1, ImIiQ{n) 
jumps to a finite value, and ReI{Q{Q) diverges logarith- 
mically H]. This behavior is shown in Fig. ga. 

Substituting S(|a;| ~ A) into (|o|) and performing the 
same calculations as before, we obtain near Ci ^1 and to 
first order in 7a 



/mnQ(ri) = —(1 + - arcsin ^ E 

2 TT I 2 



--) 



"h + l 



i?enQ(r!) = A log^ -^[a^] 

7a 



(11) 



where, as before, Qq = {^ — 1)/(7a)i and in the limits 
of small and large a, ^(a) behaves as ^(a ^ 1) = a^/2, 
and ^(|a| -^ oo) ~ log|Q;|. 

We see, similar to what we found for the Raman inten- 
sity and the DOS, the inclusion of the fermionic damping 
eliminates the singularities in the spin polarization oper- 
ator: ImllQ changes continuously through 57 = 1, and 
Re Uq is peaked but does not diverge at fj = 1 (to first 
order in 7a, the peak does not shift to higher frequen- 
cies). This behavior is shown in Fig. gb. 

At small frequencies, the expansion in f7 in (fKl) yields 



iiQin) = -ri^ + 2if„ \n\. 



(12) 



Again, similar to the result for R{i^), the inclusion of a fi- 
nite fermionic damping yields a nonzero ImIlQ{Q) down 
to the lowest frequencies. This result implies that in the 
presence of impurity scattering, the resonance peak in 
SQ{n) has a finite width. This effect may account for the 
width of the resonance neutron peak near optimal dop- 
ing. We, however, do not believe that fermionic damping 
is responsible for the increase of the peak width with un- 
derdoping - this last effect is likely to be caused by the 
frequency dependence of A associated with the pseudo- 
gap effects, which we do not consider here [ p3| . 

To summarize, we have considered in this paper a sim- 
ple model form of the electronic damping and analysed 
how it affects the forms of the Raman intensity, the DOS, 
and the spin polarization operator at the antiferromag- 
netic momentum. We found that a finite damping elimi- 
nates artificial divergencies found in a Fermi-gas consid- 
eration. Still, however, without final state interaction, 
the peak in R(fl) occurs at or beyond twice the peak 
frequency for the DOS, in contradiction with the exper- 
imental observations. This negative result implies that 



fermionic damping alone cannot account for the data, 
and supports the explanation of the downturn shift of 
the Raman peak with underdoping in terms of a mid-gap 
pseuo-resonance mode in the Raman intensity |15[ . We 
also found that fermionic damping gives rise to a finite 
width of the resonance neutron peak. 
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